Unified understanding of the valence transition in the rare-earth monochalcogenides 

under pressure 
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Valence instability is a key ingredient of the unusual properties of / electron materials, yet a 
clear understanding is lacking as it involves a complex interplay between / electrons and conduc- 
tion states. Here we propose a unified picture of pressure-induced valence transition in Sm and 
Yb monochalcogenides, considered as model system for mixed valent 4/-electron materials. Using 
high-resolution x-ray absorption spectroscopy, we show that the valence transition is driven by the 
promotion of a 4/ electron specifically into the lowest unoccupied (LU) 5d ti g band. We demonstrate 
with a promotional model that the nature of the transition at low pressures is intimately related 
to the density of states of the LU band, while at high pressures it is governed by the hybridization 
strength. These results set a new standard for the generic understanding of valence fluctuations in 
/-electron materials. 



PACS numbers: 71.20.Eh, 62.50.-p, 78.70.Ck, 75.20.Hr 

Valence fluctuations play an essential role in some of 
the /-electron systems most exciting behaviors, such as 
quantum criticality and unconventional superconductiv- 
ity [1,2]. But their understanding has proven challenging 
to capture in a unified theory because they arise from 
subtle many-body interactions between the / electrons 
and the conduction states [3, 4]. Pressure on the other 
hand is an efficient way to act on /-electron interactions 
and localization and thus can serve as a window onto 
the /-electron physics. Notably, the pressure-induced va- 
lence transition of the lanthanide monochalcogenides is 
a paradigmatic illustration of / derealization phenom- 
ena and their tremendous diversity. For instance, SmS 
undergoes a first-order transition to near-trivalency co- 
inciding with the onset of magnetic ordering [5] , whereas 
the transition of YbS into an intcrmcdiatc-valcnt cor- 
related metal is extremely sluggish [6]. These diverse 
behaviors, while establishing a severe test bed for the- 
oretical understanding of /-electron systems, have been 
overlooked for the past decades. Here, we address them 
in the light of a direct measurement of the electronic 
structure of Sm and Yb monochalcogenides (SmS, SmSe, 
SmTe, YbS, YbSe) under pressure performed using high- 
resolution x-ray absorption spectroscopy in the partial 
fluorescence yield (PFY-XAS) mode. Our data reveal 
that the valence discontinuity and concomitant closing of 



the gap under pressure are caused by the specific filling 
of the lowest unoccupied (LU) 5d ti g band by a 4/ elec- 
tron, offering a unified picture of electron derealization 
and intermediate valency. Using a promotional model, we 
show that both the steepness and the amplitude of the 
valence transition of the Sm compounds increase with the 
density of states (DOS) of this LU band. When exceed- 
ing a critical DOS, the transition becomes first order for 
SmS. The model fails to describe the prolonged transi- 
tion of the Yb compounds, which suggests that at high 
pressures the valence transition is impeded by enhanced 
hybridization. Ultimately, we suggest that our analysis 
can serve as a prototype for understanding the general 
mechanism of derealization in / electron materials. 

The PFY-XAS measurements were carried out at room 
temperature at the Taiwan beamline BL12XU, SPring-8. 
Details of the experimental setup have been published 
elsewhere [7]. For each compound, a fragment of single 
crystal of ~100 in size was loaded into the sample 
chamber of a Be gasket with silicone oil used as pressure 
transmitting medium. Pressure was achieved using a di- 
amond anvil cell. The pressure dependence of the Sm L3 
edge measured on the Sm monochalcogenides by PFY- 
XAS is shown in the top panels of Fig. 1 (a) . The spectra 
consist of two main peaks corresponding to the dipolar 
transition 2p — > 5d for the divalent and trivalent states, 
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FIG. 1. (Color online) (a) Pressure dependence of the Sm- 
L 3 PFY-XAS spectra for SmS, SmSe, SmTe (top row), and 
examples of fit performed on spectra both before (center row) 
and after (bottom row) collapse. Three components A, B, C, 
are used to fit the 2+ feature. Arrows indicate the direction 
of increasing pressure, (b) Sm 5d, t2 g , e g DOS calculations 
for SmS, SmSe, SmTe compared with the corresponding PFY- 
XAS spectra. As indicated by the solid lines, the energy shift 
between the peaks A, B, and C decreases across the series 
SmS, SmSe, SmTe. (c) Band structure of YbS. The low- 
energy, isolated part of the band A, A', is highlighted. 



respectively decreasing and increasing under pressure. A 
closer look at the leading edge of the divalent peak re- 
veals a shoulder located around 2.3 (SmS), 2.0 (SmSe), 
and 1.6 (SmTe) eV below the peak maximum (open cir- 
cle in Fig. 1(a)). A fit of the divalent component actually 
requires the use of three Gaussian functions in order to 
reproduce the low-energy shoulder, the main peak, and a 
high-energy shoulder, hereafter referred to as A, B, and 
C (cf. Fig. 1). 

In order to understand the origin of these spectral fea- 
tures, we have calculated the electronic structure of the 
Sm and Yb monochalcogenides using the density func- 
tional theory within the generalized gradient approxima- 
tion (GGA). The DOS obtained for the Sm 5d band and 
the crystal-field split t 2g and e g subbands are compared 
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FIG. 2. (Color online) (a) Pressure dependence of the va- 
lence. Yellow and cyan filled markers respectively correspond 
to before collapse (BC) and after collapse (AC), (b) Enlarged 
divalent region of three spectra per compound normalized in 
intensity to the maximum of B, measured at ambient pres- 
sure (AP), BC, and AC. The arrows illustrate the low-energy 
shift of the shoulder A between AP and BC, and its collapse 
between BC and AC. 

with the experimental spectra in Fig. 1(b) for the Sm 
compounds. We note that the Coulomb interaction is not 
taken into account and therefore the calculated bands are 
metallic. The structures A and C are found to be respec- 
tively of pure t 2g and e g character, while B is a mix of 
both. The energy splitting between the three structures 
and the relative intensity of A and C increase along the 
series SmTe < SmSe < SmS, which reflects the increase 
of the crystal field as the lattice constant decreases. We 
can therefore conclusively rule out the previous assign- 
ment of A to the quadrupolar pre-edge [5, 8]. Also, the 
band structures calculated for YbS in Fig. 1(c) and for 
the Sm compounds in ref. [9] show that the low-energy 
part of A, which we call A', corresponds to an isolated t 2g 
band sandwiched between the Fermi level and the rest of 
the conduction band, making it the LU band. 

The pressure dependence of the valence v as derived 
from the fits of the PFY-XAS spectra is shown for the 
five monochalcogenides in Fig. 2(a). In previous estima- 
tions of v, A was incorrectly assigned to a quadrupolar- 
like feature and therefore not contributing to the divalent 
weight, which resulted in an overestimation of v by ~ 0.1 
for SmS and 0.3 for YbS [5, 8]. All five compounds dis- 
play a steady increase in v at their insulator-to-metal 
transition pressure, coinciding with the pressure ranges 
reported for their volume collapse [10]. The abrupt in- 
crease in v around 15 GPa for SmTe is due to a structural 
transition. The onset pressure of the collapse gradually 
increases as a function of both the ligand (S<Se<Te) and 
more starkly the rare-earth (Sm^Yb). This reflects the 
difference in the gap values at ambient pressure, e gap — 
0.20 eV for SmS, 0.45 eV for SmSe, 0.65 eV for SmTe, 
1.4 eV for YbS, and 1.75 eV for YbSe [11-14], and the 
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idea that the wider the gap, the larger the pressure to 
close it. 

We now address the large variety of slopes and ampli- 
tudes of the transitions, and its intimate connection with 
the LU band. We start by showing in Fig. 2(b) the en- 
larged divalent region of three PFY-XAS spectra selected 
for four compounds, SmS, SmSc, YbS, YbSe, normalized 
in intensity to the maximum of B. These spectra were col- 
lected at ambient pressure (AP), just before the collapse 
(BC), and after the collapse (AC). The BC and AC pres- 
sures are respectively indicated by yellow and cyan filled 
markers in Fig. 2(a). The changes brought on by pres- 
sure on the lineshape of A are significant. First, judging 
from the shape of the shaded area between the AP and 
BC spectra, A stretches towards low energies between AP 
and BC, and the extent of the stretch increases with the 
BC pressure P BC , i.e. SmS<SmSc<YbS<YbSe. Second, 
the collapse coincides with a decrease of A, implying that 
A gets filled during the electronic transition [15]. 

In Fig. 3(b), the pressure dependence of key parame- 
ters derived from the fits of the PFY-XAS spectra gives a 
more precise view of the transition mechanism. For A, we 
consider the lower-energy portion of the peak which does 
not overlap with B, A', and the higher-energy portion 
which does, A" (cf. Fig. 3(a)), in order to distinguish 
the filling of the low-energy part of the LU band truly 
isolated from the rest of the conduction band (A') from 
that of the higher-energy, hybridized part of the LU band 
(A"). The sum Ia+Ib+Ic is scaled to the total number 
of 5d holes, i.e. between 10 for v — 2 and 9 for v = 3. 
The pressure axis is rescaled by Pec indicated by the 
vertical dashed lines. The error on the intensities and 
width is estimated to be under ±5% [15]. Confirming 
the trend observed in Fig. 2(b), a close correspondence 
is observed between the increase of v and the decrease of 
A. More strikingly for Sm, the pressure ranges of the suc- 
cessive decline of Ia' and Ia" correlate respectively with 
the valence jump and the start of the slow increase (cf. 
blue circles and black crosses in Fig. 3(b)). This demon- 
strates that the volume collapse transition results from 
the dumping of / electrons into A', while the transition 
notably slows down as soon as / electrons are promoted 
into A" , and slows even more when they start to fill the 
B band beyond AC. For YbS and YbSe the whole tran- 
sition is considerably slower, as A' is not yet filled up at 
30 GPa. The increase in Ib+c observed across the transi- 
tion for all five compounds occurs as a result of B signifi- 
cantly broadening, causing portions of A and C to merge 
with B [15]. The amplitude of the valence discontinu- 
ity is therefore determined by the remaining intensity of 
A. Also noteworthy is the finding that the monochalco- 
genide gaps are closed by a combination of both band 
broadening and shift under pressure, as seen from the 
Ea and W\ panels. 

For a systematic understanding, we have developed a 
simple model of the promotion of an electron from a 4/ 



level of width zero into a square 5d band (correspond- 
ing to the LU band) of width W and occupancy z, as 
illustrated in Fig. 3(d). W is taken as 1 eV, which is 
an approximation of Wa at BC (cf. Fig. 3(b)), and z 
is taken as vac-2- We consider that Ufd, the Coulomb 
interaction between the 4/ and the conduction electrons, 
is effective. For Ufd < |^ the valence transition is con- 
tinuous, and its pressure range, defined by the values of 
Pbc and Pac, depends on the value of Ufd- In this case, 
the energy of the Fermi level Ef is given by [15]: 



£f = 



l-2t/ /d (#) 



(1) 



where Sf is the energy of the 4/ level and E is the center 
of the square 5c? band. The valence is then expressed as: 
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(3) 



with c = e g ap /Pbc- 

For Ufd> |j , at a critical value of the pressure there is 
a first-order transition and the valence changes discontin- 
uously from 2 to 2+z. The critical pressure is determined 
based on the condition that the total energy of the sys- 
tem under pressure is equivalent to the value at ambient 
pressure. 

The transition for the Sm compounds is best repro- 
duced for the parameters indicated in Fig. 3(c). Remark- 
ably enough, the main features of the valence transition 
are properly captured for the three Sm compounds, a 
higher DOS of the LU band (^?) resulting in a steeper 
transition, ultimately first order for SmS, as shown in 
the top panel of Fig. 3(b). This lends support to the 
relevance of the promotional model, which has been an 
open question since the work of Ramirez and Falicov on 
Ce [16]. Because our model describes the linear filling of 
one band, it approximates a constant slope through the 
transition between BC and AC, and does not account for 
the filling of A" and B beyond AC. For YbS and YbSe, 
our model yields unphysical negative values of Ufd- For 
YbS for instance, using in Eqs. (2) and (3) the experi- 
mental values P B c = 8.7 GPa, P AC = 18.0 GPa, z = 0.4, 
and e gap = 1.4 cV, together with W = 1.0 eV, we ob- 
tain Ufd=-0.Q2 eV. A very plausible explanation is that 
hybridization effects between the 4/ electron and the 5d 
band become stronger at the higher transition pressures 
of the Yb compounds, so that A' no longer retains its 
isolated character. 

This analysis thus delineates two broad scenarios for 
the electronic transition of the monochalcogenidcs. In 
the first, the LU band is isolated from the rest of the 
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FIG. 3. (Color online) (a) Illustration of the decomposition of A into A' (light grey) and A" (dark grey), (b) Pressure 
dependence of v estimated experimentally (markers) and with the model (dashed line) (top row), and of the intensity of the 
peaks A', A", B+C (respectively second and third row from top). The use of different markers highlights the correspondence 
between the v transition and the filling of A', A", and B+C for the Sm compounds. Pressure dependence of the width and 
energy of the peak A (respectively second and first row from bottom). The vertical dashed lines indicate BC pressure, Pbc- 
(c) Parameters of the promotional model, (d) Schematic illustration of the promotional model. 



conduction band, allowing for a steep transition, whose 
slope and amplitude are proportional to the DOS of the 
LU band. Hybridization thrives in the second scenario, 
which leads to a sluggish transition, for which a Kondo 
description may be appropriate. We note that because 
^7 decreases and hybridization increases with increas- 
ing pressure, slower transitions may be more often found 
at higher pressures. This approximation holds for other 
rare-earth monochalcogenides, the transition being first 
order for TmSeo.45Teo.55 at 1 GPa [17], second order but 
steep for TmTe at 2 GPa[18], and sluggish between 16 
and 20 GPa for EuSe [19]. Our interpretation may also 
be generalized to other divalent 4/ materials. EuO, for 
instance has an A shoulder in its L3 spectrum, and a 
sluggish transition was recently reported above 14 GPa 
[20]. On the other hand divalent 4/ materials without 
an A shoulder show a sluggish transition which starts as 
soon as above 1 atm, see for instance Pr [21], Eu [22], Yb 
[23], YbAl 2 [24], YbGa1.15Sio.8S [7], or Cd 5 . 7 Yb [25]. We 
believe this is because / electrons are directly promoted 
into the main 5d conduction band. 

We have presented a simple framework which allows 
for a unified description of the diverse pressure-induced 
electronic phase transitions of Sm and Yb monochalco- 
genides, with implications for the 4/ delocalization phe- 
nomenon in general. We have shown that the transition 
corresponds to the filling of the LU band. This band is 
at the heart of the diverse electronic transitions of the 
monochalcogenides, as the slope and amplitude of the 
transition depend on both its DOS and degree of hy- 
bridization. Our analysis can also seemingly describe the 
delocalization in a variety of other compressed 4/ elec- 



tron materials. It may also be extended to the 5/ elec- 
tron systems, which have a broad 6d unoccupied band 
and usually show protracted electronic transitions under 
pressure [26]. 
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DATA ANALYSIS 

The pressure dependence of the Yb L3 edge measured 
on YbS and YbSe by PFY-XAS is shown in the top panels 
of Fig. SI. Examples of fits for both compounds are 
shown for the lowest and highest measured pressures in 
the center and bottom panels, respectively. 

The decrease of the shoulder A through the pressure- 
induced volume collapse, shown in Fig. 2(b), is further 
illustrated in Fig. S2. In this figure we have superim- 
posed representative spectra measured at four successive 
pressures across the collapse with the spectrum measured 
before the collapse for SmS and SmSe. All spectra are 
normalized in intensity to the maximum of the B feature. 
The shoulder A is seen to gradually decrease across the 
collapse. In Fig. S3 we show the effect on the fit of the 
±5% error which we consider for 1^ and W^. 

Fig. S4 shows the increase of Wb across the collapse. 
This broadening causes an increase of the overlap of B 
with A and C. We attribute the concomitant decrease of 
A and C illustrated in Fig. 3(a) for SmS to this enhanced 
mixing, which results in the merging of part of A and C 
with B. 
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FIG. SI. (Color online) Pressure dependence of the YD-L3 
PFY-XAS spectra measured on YbS and YbSe (top row), 
and examples of fit performed on the spectra measured at the 
lowest (central row) and highest pressure (bottom row) for 
YbS (left column) and YbSe (right column). Three peaks, A, 
B, C, are used to fit the 2+ feature. In the top-row panels 
the arrows indicate the direction of increasing pressure. 
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FIG. S2. Evolution of the feature A across the pressure- 
induced valence transition in SmS (top panel) and SmSe (bot- 
tom panel). For each compound, the S111-L3 PFY-XAS spec- 
trum measured just before collapse (BC, dashed line and cir- 
cles) is overlapped with four spectra measured during collapse 
(DC, solid line). All spectra are normalized in intensity to the 
maximum of the B feature. 



MODEL 

We first explain how the equation (1) in the article is 
derived. If we assume that n 4/ electrons per Sm are 
promoted to the 5d band in the system sketched in Fig. 
3(d), the value of n, and therefore of the Fermi level 
£p, is determined for the stable state by minimizing the 
total energy change E^and + £4/ + E int with respect to 
n (or ep): 

-, (Eband + Eif + E; mt ) = 0, (1) 

CLEp 

where Eb anc i and E^f are respectively the kinetic energy 
changes of the electrons in the 5d band and of the 4/ 
electrons per Sm, and E int the Coulomb interaction 
energy between the 5c? band electrons and the 4/ holes 
per Sm: 
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with W and z being respectively the width and the occu- 
pancy in the post-collapse state of the 5d band, and U fd 
the Coulomb interaction between the 4/ and the conduc- 
tion electrons. 

Next, the expression for the valence in the equation 
(2) of the article is obtained by: 

v = 2 + n, 

with ep and n being connected by the following relation: 
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FIG. S3. (Color online) Illustration of the uncertainty on the 
fit around the region of the shoulder A for SmS (left panel) 
and YbS (right panel). Compared with the PFY-XAS spec- 
trum (circles) are the fit (solid line), the fit where the as- 
derived value of I a was decreased (dashed line and crosses) 
or increased (dashed line and diamonds) by 5%, and also the 
fit where the as-derived value of Wa was decreased (dashed 
line and right-pointing triangles) or increased (dashed line 
and up-pointing triangles) by 5%. 
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FIG. S4. Pressure dependence of the full width at half max- 
imum of the peak B, Wb, derived from the fits for all five 
compounds. 



